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We present an approximation scheme for the calculation of the principal excitation energies 
and transition moments of finite many-body systems. The scheme is derived from a first order 
approximation to the self energy of a recently proposed extended particle-hole Green's function. 
A hermitian eigenvalue problem is encountered of the same size as the well-known random phase 
approximation (RPA). We find that it yields a size consistent description of the excitation properties 
and removes an inconsistent treatment of the ground state correlation by the RPA. By presenting 
a hermitian eigenvalue problem the new scheme avoids the instabilities of the RPA and should 
be well suited for large scale numerical calculations. These and additional properties of the new 
approximation scheme are illuminated by a very simple exactly solvable model. 
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r-| ■ INTRODUCTION 
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• Many-body Green's function theory has provided several standard approximation schemes for the calculation of 
excitation properties in atoms, molecules, and atomic nuclei. Setting out from a convenient single-particle description 



. like the Hartree-Fock approximation, these schemes usually lead to a matrix eigenvalue problem. This supplies 



approximations for the excitation energies and transition operator matrix elements of the system. The starting point 
is usually the well-known polarisation propagator fjj , a fundamental two-particle Green's function of nonrelativistic 
many-body theory. It is given by a sum of two terms that are related by symmetry (see e. g. Ref. 0]). 
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> ■ 

The so-called particle-hole part n^ r , s , (uj) already contains all the physically relevant information exhibiting single 
poles in the energy variable ui at the exact excitation energies of the system. Its name originates in the single particle 
picture which is usually taken as the zeroth order in a perturbation theoretical treatment. Zeroth order contributions 
to W h only arise if the indices r and r' are particle indices and s and s' hole indices. We speak of a particle index if it 
■ refers to a virtual single-particle state, i. e. a state which is not occupied in a Slater determinant ground state whereas 
^\ \ a hole index refers to an occupied single-particle state. The particle-hole part thus primarily describes excitations 
""^ ■ which, in a single-particle picture, may be understood as lifting one fermion from an occupied orbital to a virtual 
one. In a correlated system, however, there are no fully occupied or fully virtual single-particle states and thus the 
' ^ . exact particle-hole part contributes also for index pairs (r, s) other than particle-hole index pairs. 

The various approximation schemes may be classified by the order in which the correlation is taken into account. 
Another criterion for classification is which parts of the polarisation propagator are included. 
Ph' A simple first order approximation scheme based on W h is the so-called Tamm-Dancoff Approximation (TDA) 
jl],^. It can be seen as the first order approximation to the inverse matrix of the particle-hole part. In a quantum 
chemical language the TDA may be understood as treating the excited states on the single-excitation Configuration 
Interaction (SCI) level in comparison to an uncorrelated ground state. The Algebraic Diagrammatic Construction 
(ADC) approach represents a family of systematic higher order approximation schemes to the particle-hole part. 

Another group of approximations results from including the second part of the polarisation propagator IL,^ , s /(u>). 
This part is called hole-particle part because its zeroth order vanishes unless r and r' are hole indices and s and s' 
are particle indices. In contrast to the particle-hole part its poles appear at the negative excitation energies. The 
fundamental first order approximation scheme which treats both parts of the polarisation propagator on an equal 
footing is the so-called Random Phase Approximation (RPA) Among the many different ways of deriving 

this approximation we want to focus on the role of the RPA as the first order approximation to the inverse matrix 
of the polarisation propagator (Q). Naturally this matrix now comprises the set of particle-hole as well as hole- 
particle configurations and therefore has twice the dimension of the TDA matrix. Among the higher order schemes 
that treat both parts of the polarisation propagator we want to mention the Second Order Polarisation Propagator 
Approximation (SOPPA) [0|,§ and the Equation of Motion (EOM) method §,|l§. 
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In the Configuration Interaction picture the RPA may be seen to include ground state correlation in addition 
to the correlation of the excited states already accounted for in the TDA The RPA does so, however, in a 

non-variational manner and thus it is not obvious that the RPA yields improved results compared to the TDA. It 
has rather been observed that in certain cases the RPA excitation energies are worse than those of the TDA, also 
in systems where correlation of the ground state is of special importance (for exemplary numerical comparisons see 
e. g. Jl2|,[l3|] ) . In fact we will see later that the influence of ground state correlation onto the RPA excitation energies 
has to be regarded inconsistent with Rayleigh-Schrodinger perturbation theory (see also Ref. Q). We will further 
study a simple model system for which the RPA gives much poorer results than the TDA. For this example a new 
approximation method, called First Order Static Excitation Potential (FOSEP), yields the exact solution while posing 
a matrix eigenvalue problem of the same size as the RPA. Also a perturbation theoretical analysis shows that FOSEP 
consistently includes ground state correlation. 

The FOSEP method is a first order approximation scheme that is derived from an extended particle-hole Green's 
function Ji4| where additionally to the particle-hole and the hole-particle parts other propagator terms are present: 
The polarisation propagator has been augmented by combinations of single-particle propagators that give rise to 
additional zeroth order contributions in the particle-particle and the hole-hole index space. The introduced terms 
are chosen such that the extended Green's function now satisfies a Dyson equation which in turn defines a well- 
behaved (particle-hole) self energy in analogy to the fundamental single-particle Green's function ||l|, |l5| , |l6f . Earlier 
we have applied the general procedure to the particle-particle propagator, another well-known two-particle Green's 
function. We could show that the self-energy of a suitably chosen extended two-particle Green's function serves as an 
optical potential for elastic scattering of two-particle projectiles 0. The particle-hole self energy, in turn, may be 
understood as a sort of "optical potential for particle-hole excitations" , some general aspects of which are discussed 
in Ref. 0. 

This paper is organised as follows: After a brief review of the relevant construction principles of the extended 
particle-hole Green's functions we will define the FOSEP approximation and discuss the structure of the corresponding 
matrix eigenvalue problem in Sec. |. In the next section the properties of the FOSEP approximation for excitation 
energies and transition moments are investigated on a formal level. Fi rst th e similarity to the RPA equations is 
pointed out in order to compare the properties of the two approximations (II A). Both schemes share the fundamental 
properties of size cons isten cy ( II B ) and the invariance under unitary transformations of the occupied or v irtual 
Hartree-Fock orbitals ( jfl C| ) . The differences between the two approximations will become apparent in Sec. [IE when 
we will perform a perturbation theoretical analysis of the excitation energies up to second order. This analysis shows 
that the FOSEP approximation includes part of the ground state correlation in a consistent way while the RPA 
proves inconsistent with perturbation theory in this respect. It follows an analysis of the approximation for transition 
moments which is found to be consistent in first order. The last paragraph of Sec. [il] deals with the equivalence 
of length and velocity form of the transition moments of the dipole operator which are important for the so-called 
oscillator strengths. In Sec. Ill the FOSEP method is compared with the RPA and TDA in application to a simple 
exactly solvable model system. 



I. THE FOSEP APPROXIMATION 



The extended particle-hole Green's function is one species in a family of two-particle propagators that fulfil Dyson's 
equation. The general theory has been developed in Ref. |l4|,[l7]]. Here we will only outline the main ideas that are 
relevant in the present context. The desired physical information like excitation energies and transition moments is 
contained in the poles and residues of the extended Green's function. Calculating these is equivalent to solving the 
eigenvalue problem related to a generalised excitation energy operator Ti. which lives in an extended Hilbert space Y. 
Throughout this paper we assume a discrete eigenvalue spectrum for the relevant Hamiltonians since we are interested 
in finite basis set approximations. The generalisation to continuous spectra, however, poses no particular problem. 

A basis of the extended Hilbert space Y is given by the set of states { \Qi)} which are chosen "orthonormal" with 
respect to an indefinite metric ft: 

(Qi\(i\Qj)=±Sij (2) 

A model space is spanned within the full space by a subset of basis states {|Y~ rs )}. This model space is supposed to 
allow the description of those excitations that are predominantly of particle-hole type. Direct reference to 'occupied' 
and 'virtual' single-particle states is avoided by allowing the indices r and s to range over the full set of single-particle 
indices each. The states \Y rs ) are chosen to include (ground state) correlation being 'correlated excited states' in the 
sense of Refs. [ |l8||l9f| . On the other hand they are constructed manifestly 'orthonormal', i. e. satisfying 

(Y rs \iJ,\Y r > s >) = 5 rr i8 ss ' (3) 
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exactly and in each order of perturbation theory. It is a special property that all states in the primary subset {l^rs)} 
have positive norm. The construction of states with the described properties presents the crucial step in developing 
the theory of extended particle-hole Green's functions. Explicit expressions for the particular choice used in this 
paper can be found in Rcf. |l4| together with a thorough discussion of the construction principles and the remaining 
freedom of choice. 

The basis { \Qi}} defines a matrix representation H of the generalised excitation energy operator Ti . The subdivision 
of the basis into the basis of the model space {l^s}! an d the complementary part superposes a block structure onto 
this matrix: 

The index a refers to the model space and b to its complement. The primary block TL of this matrix is given by 



H 



=a a 



{Y rs \(M\Y r , s ,). (5) 



Explicit expressions for this matrix may be found in the appendix and are derived in Ref. p4j where also a physical 
interpretation as a static particle hole scattering potential in the case of Coulomb interacting particles is given. In 
the following we will investigate the first order approximation to this matrix within the framework of many-body 
perturbation theory. 

In order to apply perturbation theory the many-body Hamiltonian H has to be split into two parts H Q and Hi as 
usual. The choice of the one-particle operator 

H = ^^eia\a,i (6) 

i 

as characterised by the diagonalising single-particle basis { \tpi}} and single-particle energies {si} defines the zeroth 
order of perturbation theory. The residual interaction Hi = v + V contains contributions of a one-particle operator 
v and a two-body interaction V: 

V = i ^2 a\ataia k (8) 

In particular we are interested in the M0ller-Plesset partitioning of the Hamiltonian where the zeroth order Hamil- 
tonian H Q is defined by the Hartree-Fock approximation. In the case of a nondegenerate ground state the matrix 
elements Vij of the one-particle part of the interaction are then given by 

vf/ = -J2 n ^m- (9) 

k 

Here V rs t r ' s '] — V rsr ' s i — V rss i r i denotes the anti-symmetrised matrix element of the two-body interaction and n r is 
the occupation number of the orbital \tp r ) in the zeroth order ground state Slater determinant |$ ). 

The zeroth order of the primary block H from Eq. (0) yields the matrix e of zeroth order excitation energies: 



H 



(o) 

, , = Mrs rV = ( £r ~ £fi ) S rr' S ss> (10) 



This reflects the fact that in our ansatz the first index in the pair rs has to be understood as marking the orbital 
(position) into which a particle is created and the second index as marking the orbital where a particle is destroyed 
(or a hole created). 

The main motivation for developing the theory of the extended particle-hole Green's function was that it fulfils a 
Dyson equation and therefore possesses a particle-hole self energy S_(u>) [|l^,^7[ . This self energy may be seen to result 
from a partitioning of the eigenvalue problem associated to the matrix "H with respect to the primary block H aa - The 
energy independent ("static") part of the self energy 5(oo) is defined by the primary block minus its zeroth order: 

S(oo)=H -e (11) 
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The significance of this part is to describe the the influence of correlation to particle-hole excitations that remains in 
the high energy limit, i. e. when the "target particles have no time to rearrange upon the influence of the particle hole 
excitation" . The lowest order contributions to this matrix are of first order. The energy dependent ( "dynamic" ) part 
of the self energy S_(ui) — S_(oo) takes account of the remaining blocks of the matrix H and starts in second order. 
The first order contributions to the (static) self energy are given by 



[£M] ,, 

L— \ / j rs r ' s 



(i) r id) 



H 



V rr '8 ss > V s / S $ rr t 

+V rs '[ sr i] (n r n s - n r n s )(n r >n s i - n r >n s ,) 

+S SS > ^ n kVrk[r'k] - 5rr> ^ n k V s'k[sk] (12) 
k k 



for a general Hamiltonian where we have introduced the notation n r = 1 — n r . Expression (|12|) may be readily derived 
from the general expression for the primary block H aa which can be found in the appendix. In the M0ller-Plesset 
case ([)]), the first order simplifies even further. In particular it is this approximation that we will refer to as First 
Order Static Excitation Potential (FOSEP): 



= -I rs,r's' 



= V rs ,[ sr ,} (n r n s - n r n s )(n r ,n s > - n r m s i) (13) 

Now the following hermitian eigenvalue problem remains to be solved: 

{g + g FOSEP ) X= ujX (14) 

The (physical) eigenvalues u> provide approximations to the excitation energies of the system and the corresponding 
eigenvectors X_ may be used to calculate transition operator matrix elements (transition moments). The transition 
moments corresponding to the dipole operator define the so-called oscillator strengths which are of great importance 
for photo-absorption and emission processes j|] . The FOSEP approximation for transition operator matrix elements 
reads @ 

(tf |f|^r osBP = ^T«X£., (15) 

ij 

where T JJ are the matrix elements of the (one-particle) transition operator T and X£ are the components of the 
eigenvector associated to an excitation into the state |'I' AI ). 

Due to extensions included in the definition of the primary states |^ s ); however, not all of the eigenvalues and 
eigenvectors of Eq. (|lj) correspond to "physical" excitations. More insight may be gained by looking at the particular 
block structure of this eigenvalue problem which is revealed when splitting the set of index pairs according to whether 
the indices relate to occupied (hole) or unoccupied (virtual or particle) orbitals. In zeroth order only the diagonal 
matrix e is present and the distinction between "physical" and "unphysical" excitations is obvious, since the ground 
state is approximated by a Slater determinant. In this case excitations of a single particle are possible only from an 
occupied into a virtual orbital. Thus only the ph-ph block is "physical" . In the first order (FOSEP) approximation 
the secular matrix has the following structure: 
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e pp 



(16) 



The asterisk (*) denotes complex conjugation for the matrix elements and the tilde (~) denotes a simultaneous trans- 
position of the two index pairs that specify a matrix element which simply means a renumbering of the rows and 
columns of that matrix: 

r A 

First of all we notice that the FOSEP self energy S_ FOSEP does not contribute at all for pairs of orbitals that are 
both occupied (hh) or both virtual (pp) and only the zcroth order matrix e remains. Thus the hh-hh and pp-pp blocks 
decouple from the rest of the matrix and the eigenvalue problem ([fl]) for these blocks becomes trivial, simply yielding 
the Hartree-Fock orbital energy differences. These blocks are obviously not correlated in the first order treatment. 
We want to mention that this decoupling of the hh-hh and pp-pp blocks is special to the M0ller-Plesset partitioning 
of the Hamiltonian and to the particular choice of the extended states \Y rs ) considered in this paper. It does not 
appear for other choices discussed in Ref. Q . The decoupling leads to a considerable reduction of numerical effort 
and therefore justifies the present choice. Many of the properties discussed in the present paper, however, generalise 
also to first order approximations based upon other choices for the primary extended states |Y rs )- 

The ph block of the e matrix e ph contains those energies that relate to a simple particle-hole excitation in a zeroth 
order picture. The contribution of 



ph,p' h' 



Vph'[hp'] 



(h, h' occupied, p, p' virtual orbitals) 



(17) 



in the ph-ph block describes the interaction of the uncorrelated ground state Slater determinant with a singly excited 
configuration. In fact, diagonalising the ph-ph block e ph + V_ on its own results in the well-known Tamm-Dancoff 
Approximation (TDA) Q. The coupling W to the hp-hp block can be understood as taking into account ground state 
correlation as will be explained later with the help of perturbation theoretical arguments. A similar coupling also 
appears in the Random Phase Approximation (RPA) Q. The relation of our approach to the RPA will be discussed 
in detail below. 

The hp-hp block by itself does not seem very physical in character at all. Its zeroth order excitation energies are 
negative and result from creating a hole in a virtual orbital and a particle in an occupied orbital. Nevertheless, the 
hp-hp block couples through the matrix 



m 



ph.h'p' 



= -v. 



pp 1 [hh 1 



(/i, h' occupied, p, p' virtual orbitals) 



(18) 



with the physical ph-ph block and thus introduces a correction to the Tamm-Dancoff excitation energies. Due to 
the decoupling of the pp and hh blocks we are left with an eigenvalue problem comprising the blocks of the FOSEP 
matrix with ph and hp indices: 



e + S f 



s ph _|_ y w 



W -e ph + V 

ph and hp blocks \ — ^ 



(19) 
(20) 
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This is the eigenvalue problem that has to be solved in the FOSEP approximation scheme. In contrast to the RPA it 
is a hermitian eigenvalue problem always yielding real eigenvalues. In fact, in most cases the matrix M_ FOSEP is real 
symmetric. 

As long as the interaction remains weak enough, there is a clear distinction between "physical" eigenvalues of this 
matrix and "unphysical" ones by the sign of these energies. Even with stronger interaction the distinction may still 
be valid. This is understood easily when considering a model system where the ph-ph block has the dimension one 
and all matrix elements are real. The matrix ]\£ FOSEP then is two by two and its eigenvalues are given by 



±\J eP h2 + W 2 + V. 
Thus we get one positive eigenvalue and a negative one, provided that 

V 2 < e ph2 + W 2 . 

Note that V is the contribution of the correlation introduced in the TDA and this condition now states that it has 
to be small enough compared to the zeroth order excitation energy augmented by the additional interaction term W. 
When this condition is violated, or in general the numbers of positive and negative eigenvalues in a given symmetry 
are not the same, the eigenvectors may be necessary to distinguish between physical and unphysical contributions. 
Still the "physical" approximation may usually be defined by the upper half of the eigenvalues. 

In the remaining sections of this paper we will discuss the FOSEP approximation as defined above. It presents 
the natural first step in approximating the particle-hole self energy S_(u>) and thus the matrix "H. At this place we 
want to mention that other approximations, for example, result by augmenting the primary set of states \Y rs ). In 
particular one can obtain the RPA straightforwardly via the formalism of extended Green's functions. This is achieved 
by additionally including a subset of the basis { \Qi)} consisting of states with negative norm that are degenerate (in 
zeroth order) to the hole-particle fraction of the set {|^rs)} fb|- This augmented set of states defines an extension 
of the primary block TL in the matrix "H. In first order this extended matrix can be decoupled with the help of a 
unitary transformation into the RPA eigenvalue problem and additional unphysical blocks. Thus the RPA is included 
in the general theory as a specific approximation. Note, however, that FOSEP presents the canonical first order 
approximation in our ansatz since it is based upon the primary set of states {|i^ s )} which defines the Green's function 
and self energy matrices. 

II. PROPERTIES OF THE FOSEP APPROXIMATION 

In this chapter we will discuss some general properties of the FOSEP approximation. In order to elucidate the 
relation to the well known first order approximation schemes RPA and TDA we start with briefly reformulating the 
RPA in our notation. We go on with considering two fundamental invariances FOSEP shares with RPA and TDA 
namely size consistency and invariance with respect to unitary transformations of the single-particle basis. In order to 
show the differences between the three schemes we will carry out a perturbation theoretical analysis for the excitation 
energies as well as for the transition moments. Finally, the equivalence of length and velocity form of the dipole 
operator transition moments is discussed. 

A. Relation to the RPA 

The Random Phase Approximation (RPA) QJ^] for the calculation of excitation energies and transition operator 
matrix elements in finite Fermi systems may be derived and understood in many different ways. Traditionally the 
RPA is derived by the infinite summation of a certain type of diagrams in the Feynman-Dyson perturbation series 
of the polarisation propagator || . Equivalently it can be understood as a first order approximation to the integral 
kernel of the Bethe-Salpeter equation pofl or as a specific first order approximation in the equation of motion of the 
polarisation propagator . Now we are going to present the RPA equations in a form suitable for comparison with 
FOSEP. 

Based upon a Hartree-Fock zeroth order the RPA is defined by the following eigenvalue problem [Ml : 

(em + P^x = u>mx (21) 

The matrix s is defined like above in Eq. (|l0|). The RPA kernel R consists of the matrix elements R rs ,r's' = Ks'[sr'] 
and the metrical matrix to by 



G 



\m] , , = S rr 'S ss i (n r n s — n r n s ). (22) 
Note that in block matrix notation m can be written as 



I g 
-1 



(23) 



All these matrices are indexed by pairs of single-particle indices that are required to be either particle-hole or hole- 
particle index pairs. Thus the RPA eigenvalue problem has the same size as the FOSEP one (|l9|). The essential 
difference is the appearance of the indefinite metric m in the RPA case which renders the RPA problem a non- 
hermitian eigenvalue problem. As a consequence the RPA may become instable and lead to complex eigenvalues [§). 
The RPA kernel E is related to the FOSEP self energy £ FOSBP m the following way: 

g FOSEP =^Em (24) 

Introducing the matrix 



M RPA = em + S FOSEP 



T_ph + y 



(25) 



ph and hp blocks 

with the same nomenclature as in (|l^) and (|l8|), the RPA eigenvalue problem may be rewritten to 

M Z RPA x' = umx! (26) 

with x' = mx. Comparing the RPA (^6|, ^5j) with the FOSEP eigenvalue problem (|T|, |2(]) we see that both have 
the same size and start from the same input data while the difference lies in some minus signs. Before analysing the 
differences further we want to discuss two fundamental properties that are shared by both schemes. 



B. Size Consistency of FOSEP 



The question of size consistency of a many-body method is the question of whether the resulting approximations 
for physical quantities scale correctly with the size of the system (2^j2j| . The general question is difficult to answer 
and usually one has to resort to simple models or numerical calculations. Nevertheless, this concept becomes very 
important for applications to large or extended systems. In the context of finite systems, especially molecules, the so- 
called separate fragment model provides a useful test of correct scaling behaviour. We consider a many-body system 
consisting of two or more separate (noninteracting) subsystems (fragments). Size consistency of excitation energies 
and transition moments then means that an excitation that is local to one of the subsystems is approximated with 
the same result regardless whether the approximation scheme is applied to the full system or only to the fragment. A 
sufficient but not necessary condition for this property is that the secular equations of the approximation scheme give 
rise to independent sets of equations corresponding to local excitations on the individual subsystems. This 'a priori' 
decoupling of independent, local equations is known as the separability property p4| . 

The separability of FOSEP and RPA is proven by the following arguments: In the model of separating fragments 
the Hamiltonian of the full system is given by the sum of the Hamiltonians of the subsystems. This implies that the 
(Hartree-Fock) single-particle states <p r may be chosen local to either fragment and that the matrix elements of the 
two-body interaction Vijki vanish unless all indices ijkl relate to states belonging to the same subsystem. From the 
definition of the FOSEP matrix (]l3|, ^) and the RPA matrix ( p5|) it therefore becomes clear that both methods have 
the separability property and thus can be regarded size consistent. The same arguments apply to the TDA to which 
the FOSEP and RPA reduce in the case of vanishing coupling W_ as explained in Sec. |. For nonlocal excitations the 
excitation energy is simply given by the difference of the single-particle (Hartree-Fock) energies for all three schemes. 
This means that the level of approximation is that of Koopman's theorem which provides a consistent first order 
description. We want to mention that the separability property is by no means a matter of course for more accurate 
many-body methods like, for example, the general CI method p5[ . 



C. Unitary Transformations of Single— Particle Orbitals 



In order to separate the influence of finite truncations of the underlying single-particle basis from the systematic 
deficiencies of a given approximation scheme it is important that the approximation is invariant with respect to 
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rotations of the single-particle basis. A global invariance is self-evident only for 'exact' methods like full CI. Systematic 
truncations of the CI matrix employing single, double, or triple excitations on a given reference configuration are at 
least invariant with respect to transformations of the single particle basis that do not mix occupied and virtual orbitals. 
Such an invariance usually does not apply for perturbative propagator methods. The FOSEP method as well as the 
RPA and the TDA, however, share this invariance with the CI whereas higher order methods usually do not. For the 
SOPPA method the influences of rotations of the orbital set have been investigated numerically |26[ . 

A physical motivation for altering the single-particle functions may be drawn from the fact that the Hartree-Fock 
virtual orbitals describe additional test particles in a mean field and therefore constitute rather diffuse functions while 
the main effects of correlation show up at short range due to ineffective screening. Thus one can hope to achieve a 
better description of the influence of correlation with more localised virtual orbitals than the Hartree-Fock ones. 

Within a perturbation theoretical approach a unitary transformation within the set of virtual single-particle states 
can be realised by adding a (hermitian) single-particle potential to H of Eq. (^) that takes effect only on the virtual 
orbitals and subtracting it again fromiJi of Eq. (@). The new single-particle basis is then defined as the diagonalising 
basis for the new zeroth order Hamiltonian. Obviously the new basis is connected to the original one by a unitary 
transformation that leaves invariant the occupied single-particle states and also does not affect the Slater determinant 
|$ ) preserving the distinction between occupied and virtual single-particle states. From the definitions of the matrices 
( p"o| , |l3| ) it can be seen that such a transformation of the single-particle basis also results in a unitary transformation of 
the secular matrix in (|lj) which preserves its block structure (|l6|). Therefore also the FOSEP eigenvalue problem ( fisj , 
pfj| ) transforms without changing its eigenvalues. This argumentation can be transferred analogously to the closely 
related RPA and TDA. Summarising, we have seen that the FOSEP approximation as well as the RPA and TDA 
are invariant under unitary transformations within the set of virtual single-particle states. It is easily seen that this 
property generalises to unitary transformations of the orbital basis that do not mix occupied and virtual orbitals. 

In a much more general sense, however, the matrix H aa of Eq. (^|) which forms the primary block of the matrix 

representation 7£ of the excitation energy operator Ji is invariant under (un rest ricted) unitary transformations of the 
single particle space. Since the only unknown quantity in 7£ (see also Eq. (Al) in the appendix) is the exact ground 
state of the system, the invariance properties of an approximation to iL aa follow the chosen approximation for the 
ground state. In other words: The eigenvalues of the matrix of Eq. (Al) depend only on the chosen approximation 
for the ground state and not on the particular choice of the single-particle basis. In the special case of a system of 
particles which interact only with a one-particle potential, the primary block 7£ is even independent of the ground 
state. In this case the first order already provides the exact solution for the excitation energies and is invariant 
with respect to any unitary transformations of the single-particle basis. This property is explained in more detail in 
Ref. [H. 



D. Perturbation Theoretical Analysis of the Excitation Energies 



In order to analyse the differences between the FOSEP approximation, the RPA, and the TDA we will now perform 
a perturbation theoretical analysis of the excitation energies up to second order and compare with straightforward 
Rayleigh-Schrodinger perturbation theory following Ref. |Q. Usually the Rayleigh-Schrodinger series itself is not a 
reliable method for calculating energies of excited states but it is very helpful for analysing and comparing different 
approximation schemes. 

We assume that the Rayleigh-Schrodinger series starting from the singly excited Slater determinant 

|* Q/9 ) =aia p \<P ) (27) 

converges towards the excited state \^ a /3)- Note that here a has to be a particle index and j3 a hole index. 

An expression for the second order excitation energy can be gained by subtracting the expressions for the ground 
state energy Eq from the excited state energy E a p. Up to second order the ground state energy is given by the familiar 
expression 

E (2) = E a (l) + uf p - 2h \ (28) 
where E (l) = ($ | H |<J> ) is the first order ground state energy. The term 

TT (2p-2h) /on\ 

C/o = - > , n inj n k ni 29) 

j^j Ei+ej - e k - et 

k<l 
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denotes the second order contribution to the ground state correlation. The given approximation is known as the 
M0ller-Plesset (MP)-2 approximation and extensively used in quantum chemistry. In the Configuration Interac- 
tion language the second order term Ug 2p 2h%> may be interpreted to present interactions of the ground state Slater 
determinant |$ ) with two particle-two hole configurations. 

The second order energy of the excited state is also evaluated straightforwardly and can be found in Ref . . The 
excitation energy up to second order then reads 

AE a p(2) = E a g(2) - E (2) 

= AE af3 {l) + U%- h) + U^ 2h) + R a p (30) 

where AE a p(l) = e a — sp — V a p[ a p] is the first order excitation energy. The terms U^f^ and U^ 2p ~ 2h ^ denote second 
order contributions to the excited state's energy arising from the interaction of the configuration |$ Qj a) with (other) 
p — h and 2p — 2h configurations, respectively, and can be found in Ref. [Q]. The part R a p is the remainder of a 

partial cancellation of the second order ground state correlation term Uq 2p 2h ^ of Eq. ( p9f ) with a contribution to the 
correlation of the excited state. It can be written as the sum of three terms 

R a p = Rip + R 2 a p + R a p (31) 



where 

R *P = 



\Vaj[kl]\ 



k<l 



E a + Sj — £fc — £; 



R% = £ ~ - ^n, (32) 



2 



R aP - Zl 



ajlPl] 



2 



e a + Ej — ep — ei 

These contributions are left over from Uq 2p 2fl ' corresponding to the special cases where i = a or fc = in the sum 

of Eq. (|2^). The rest of Uq 2p is cancelled by contributions from the excited state. 

We are now in the position to compare with the second order excitation energies from the TDA, the RPA, and 
the FOSEP scheme. The approximations for the excitation energies in these schemes are found by solving the 
eigenvalue problem related to the corresponding matrix. Basic matrix perturbation expansion leads to the second 
order approximation for the eigenvalue. We find the following: 

AE^ A (2) =AE ad (l) + U%- h) (33) 
AE^(2) = AE aa (l) + U {p ; h) - R% (34) 
AE^ P (2)= AE aa (l) + U (p - h) + Rl B (35) 



All three approximation schemes are consistent in first order with the Rayleigh-Schrodinger expression (30). Therefore 
they are correctly referred to as first order schemes. The zeroth and first order contribution AE a p(l) originate in 
the diagonal matrix elements of the ph-ph block (TDA block) of the matrices Mf OSEP and M? PA of Eqs. (||, ||) 
whereas the second order terms come in by the diagonalisation procedure. Neither of the three schemes reproduce 
the second order expression ([SO]) completely. This is only achieved by more accurate and more costly schemes like 
ADC(2) or SOPPA. 

The term U^L ^ describes part of the second order correlation of the excited state as can be seen from Eqs. ( |30| ) 
and ( p8| , p9| ) . It is the only second order contribution to the TDA excitation energy and originates in the off-diagonal 
part of the TDA matrix which also constitutes the ph-ph block of the FOSEP matrix Mf OSEP . In the FOSEP and 
the RPA expressions the additional term i?^ is present which has already been identified as part of the ground state 
correlation. Comparing with the second order perturbation expansion (|30|), however, we notice that it appears with 
the wrong sign in the RPA excitation energy in contrast to FOSEP where the sign is consistent with the Rayleigh- 
Schrodinger expansion. This term comes in through the coupling of the ph with the hp block in the FOSEP matrix 







and the RPA matrix M_ RPA of Eq. ( p5| ) , respectively. We recall that the only difference between 
the FOSEP and the RPA equations are minus signs in the hp-hp blocks of the secular matrices and the RPA metric. 
These carry through to the second order expressions and prove inconsistent with Rayleigh-Schrodinger perturbation 
theory in the RPA case. 

It is interesting to note that each of the three parts of R a p of Eq. ( |32| ) is positive. Thus the RPA always lowers 
the TDA excitation energies in second order whereas the FOSEP approximation increases the energies in accordance 
with the positive sign of the full correction R a p. Taking into account that the TDA approximates the ground state by 
the variational Hartree-Fock method, it is sensible to expect an increase in the excitation energies when ground state 
correlation is additionally taken into account. Therefore we may conclude that the FOSEP scheme includes ground 
state correlation in a consistent way in contrast to the RPA which does not. 

Concerning the term U^~ 2h ^ in the full second order expression ( |30| ) which is missing in all three schemes, we 
remark that it does not carry a definite sign (as can be seen from Eq. (68c) of Ref. J2j). It may, however, lower the 
energy of the excited state and thus can possibly overcompensate for the influence of the ground state correlation. In 
this way the RPA result may be supported by accidental numerical compensation. 



of Eq. (19 



E. Perturbation Theoretical Analysis of the Transition Moments 



We now will show that transition operator matrix elements between the ground state and an excited state can, in 
contrast to TDA, be approximated by FOSEP consistently in first order. We consider the transition moment 

T a p = {^ \f\^ a p) (36) 

of the (one-particle) transition operator T for the particle hole excited state \^ a p) that was introduced in the last 
paragraph. Up to first order the perturbation expansion yields 

r a(3 (l) = f \* afi ) + ($ | f |¥<g) + ^ { o ] \T |$ Q/3 ). (37) 

Explicit expressions for these terms derive from straightforward Rayleigh-Schrodinger perturbation theory and may 
be found in Ref. §. 

In an analogous fashion to the preceding paragraph the above expression may be compared to the result of matrix 
perturbation analysis of the TDA, RPA, and FOSEP approximations for the transition moments. It is easily seen 
that the FOSEP eigenvalue problem (jl^, |l^) together with the approximation for the transition moments of Eq. ( |T5| ) 
leads to a consistent first order approximation of the transition moments. For the TDA, the RPA, and the first and 
second order ADC scheme the perturbation analysis has been carried out in Ref. Q]. It has been found that the 

TDA expression is incomplete in first order because the term (^o^l T |$ Qj a) is missing. This term is a consequence of 
first order ground state correlation which is neglected in the TDA. It is described correctly in the RPA which yields 
consistent transition moments through first order. This may be seen to justify the common statement that the RPA, 
in contrast to the TDA, considers ground state correlation. We want to mention that the transition moments are 
also described consistently through first order in the ADC(l), an approximation that is equivalent to the TDA with 
respect to the excitation energies but differs for the transition moments 0] . 



F. The Equivalence of Length and Velocity Form of the Dipole Operator Transition Moments 

Now we will focus on the transition moments of the dipole operator as a particular choice for the transition operator. 
These transition moments are related to the so-called oscillator strengths and present important parameters charac- 
terising the interaction of the many-body system with radiation. There is a hierarchy of equivalent representations 
of the exact dipole operator transition moments, starting with the so-called length and velocity form. Without loss 
of generality we consider only the ^-component Z of the dipole operator. The identity 

- Eo)(*p\ Z |tf ) = -1(^1 A (38) 

expresses the equivalence between the length and the velocity form of the dipole operator transition moments defined 
by the left and the right hand side of this equation, respectively. The z-component of the momentum operator P z is 



related to the dipole operator by H , Z 



-iP z provided that the Hamiltonian contains only local potentials like e. g. 
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for Coulomb interacting electrons in atoms or molecules. Eq. (|3q) then follows from the identity (E^—Eq)^^] Z \9 



H, Z 



I*, 



It is a very special property of the RPA to preserve this equivalence exactly, provided the underlying Hartree- 
Fock single-particle basis is complete p7[| . In fact the RPA eigenvalue problem may be derived setting out from a 
CI representation of the ground and excited state wavefunction and requiring certain 'hypervirial relations' which 
present a slight generalisation of the above mentioned equivalences [[ll) . Thus the exact fulfilment of these hypervirial 
relations may be regarded unique to the RPA while other approximation schemes show, at most, a perturbative 
equivalence. In order to make the perturbation expansion transparent, it is useful to introduce the function A(A) as 
the difference between the right and the left hand side of Eq. (pq ) . It is a function of the usual interaction strength 
parameter A in H\ = Hq + \H\. Owing to its non-local character, Hartree-Fock does not preserve the length- velocity 
equivalence. Hence, choosing Hq to be the Hartree-Fock operator, A(A) does not vanish for any A except for A = 1. 
The expansion of A (A) in powers of A does not vanish term by term; each order will in general give a nonvanishing 
contribution. The TDA approximates only the zeroth order term of A(A) correctly, giving an error in first order [i. 
e. 0(A)]. The FOSEP approximates A(A) through first order giving a second order error. So does the RPA which 
additionally has the unique property of reproducing the zero A(A = 1) = 0. For completeness we mention that the 
first order ADC j28| approximates A(A) through first order while the SOPPA |^9| as well as the second order ADC 
p8[ are consistent through second order. 

The result of equivalence through first order is another clue for the consistency of the FOSEP approximation. The 
error between length and velocity form provides a genuine test for the validity of the approximation. Depending on 
the particular aims of an approximate calculation this may seem more favourable than having an 'a priori' equivalence 
like in the RPA which hides the actual error with respect to the exact value of the transition moments. 



III. APPLICATION TO A SIMPLE MODEL 



In this section we compare the FOSEP approximation with the RPA, the TDA, and the exact solution for the 
excitation energies of a very simple model system. This model is commonly referred to as the Hubbard model for 
the Hydrogen Molecule H2 [^0|. In the model, all excitation energies can be calculated by analytic expressions in 
dependence of two parameters which mimic the effects of Coulomb interaction. 



A. Definition of the Model and Exact Solution 



The Hubbard model for the diatomic Hydrogen molecule represents each atomic site (R and R') by a single orbital 
electronic level denoted by |R) and |R'), respectively. Each level can be occupied by up to two electrons with opposite 
spin. The single-particle part of the Hamiltonian h contains a diagonal term 

(R|fc|R) = <R'|fc|R')=£ (39) 

which yields an energy £ for each electron. The off-diagonal term describes attraction by the neighbouring nucleus 
and represents an amplitude for 'tunnelling' or 'hopping' of an electron from one site to another: 

(R|fc|R') = <R'|&|R)=-t (<0) (40) 

Additionally a two-particle interaction term is present which yields a positive energy U whenever one level is occupied 
by two electrons. This term represents the intra-atomic Coulomb repulsion between two localised electrons. Note 
that all interactions are independent of electron spin. 

We consider the neutral H2 molecule, i. e. an occupation with two electrons. The solution of the Hartree-Fock 
equations yields the (molecular) orbital functions 

= -L( |R> + |R'» |«) = ^=( |R> - |R')). (41) 

while the corresponding Hartree-Fock single-particle energies read e g / u = E^ft+^U. In the Hartree-Fock ground state 
Slater determinant |$ ) = \g ] g {) the orbital \g) is doubly occupied. In order to make transparent the transition 
between the Hartree-Fock approximation and the correlated problem we introduce the additional (perturbation) 
parameter A G [0, 1] by using the Hamiltonian 
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H X = H + XHi (42) 

where the Fock operator H and the interaction Hamiltonian Hi are defined as in Sec. |J The matrix elements of the 
two-body interaction Vijki that define Hi are given by the transformation ( fill ) into the atomic orbital picture. 

There are six independent solutions for the two-electron eigenstates of the system. According to the possible 
combinations of the electrons' spins three states of singlet symmetry l^o/i) an d \S U ) and a degenerate triplet |r_ 1/0/1 ) 
are found. The triplet and the singlet \S U ) states are uncorrelated and have one electron in a u and one in a g orbital. 
The singlet state |<J/i) corresponds to a 2p — 2h excitation. It will no longer be considered because it lies outside of 
the range of RPA, TDA, and FOSEP. The wavefunction of the singlet ground state is given by 



\V }x=a x \gU D+Px\uUl) (43) 



where 



ax 



4t + Vl6< 2 + A 2 U 2 



yj\ 2 U 2 +(4t + Vl6t 2 + \HP) 2 ( 44 ) 



Pi = l-o£. 

The dependence on the perturbation parameter A indicates the influence of correlation which only appears between 
the singlet Slater determinants of g symmetry \g f g {) and \u f u I). In order to simplify the notation, we will 
drop the subscript A in the following. Note that the uncorrelated case corresponds to A = where a = 1 and (3 = 0. 
Therefore it is clear that \^f a ) is connected to the Hartree-Fock ground state |$ Q ) = \g 1 g J.). The energy eigenvalues 
of the exact states are given by 



E = 2£ + U-±U- ^jAt 2 + ^U 2 

E Su = 2S + U ' ( 45 ) 

E T = 2E + U-XU 

The excitation energies AEi are defined as usual by the difference of the excited state's energy Ei to the ground state 
energy Eq. The energy related to the excitation into the triplet, e. g., is thus given by 



A / A 2 

AE T = E T -E = --U+J4:t 2 + — U 2 (46) 
Note that the expansion of AEi into a power series in A yi elds the Rayleigh-Schrodinger series of the excitation 



energies which has been discussed in a general context in Sec. [ID . 



B. Results for TDA, RPA, and FOSEP 

We now discuss approximations to the (singlet) excitation from |^ ) into \S U ) and the (triplet) excitation into one 
of the \Ti) states. E. g. the triplet excitation into |T_i) is defined by a particle-hole excitation from the \g J.) to the 
| it t) orbital. The TDA matrix for this triplet excitation is one dimensional because the excited state is uncorrelated. 
Therefore the TDA result for the triplet excitation energy 

AE^ DA = 2t- -U (47) 

coincides with the first order of Rayleigh-Schrodinger perturbation theory. 

The FOSEP matrix M_ FOSBP as well as the RPA matrix, however, have dimension two because of the coupling with 
the corresponding hole-particle configuration. The FOSEP matrix for the triplet excitation is given by 

M FOSEP = ( 2t 7 z U 2 U \ ( 48 ) 
=t \ |f/ -2t-$U J 



Its eigenvalues are given by 



ae fo SBP (pM =-±U± JaV + ^TJ- (4!)) 
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As discussed in Sec. | we obtain physical and unphysical eigenvalues of which the latter do not carry any physically 
relevant information. The distinction between the physical and the unphysical eigenvalue is clear in the present case 

is non-positive for all choices of the parameters. Note 



because AE, 



FOSEp(p) 



is always non-negative while AE^° SEP ^ 



that the physical eigenvalue yields the exact excitation energy AE^ OSBF ^ p ' — AEt\ 

The close relation of the RPA eigenvalue problem to the FOSEP matrix was discussed in Sec. II A. The solutions 
of the RPA equations for the triplet excitation are given by 



AE, 



rpa(p/u) 



±V4i 2 -2\tU 



(50) 



Obviously the expression under the square root may become negative for certain choices of the parameters, in which 
case the RPA becomes unstable. 

A plot of the solutions of the FOSEP, the RPA, and the TDA equations as a function of the perturbation parameter A 
for a particular choice of the Hubbard parameters t and U can be found on the top of Fig. 1 . For this choice representing 
strong interatomic Coulomb repulsion, the RPA becomes instable. As discussed in Sec. II D in the framework of a 



second order analysis, the RPA lowers the TDA value whereas the FOSEP approximation yields a higher value for 
the excitation energy which is correct in the present case. We have already mentioned that here only the ground state 
|\I/ ) is correlated whereas the excited state |T_i) is not. This ground state correlation is correctly taken into account 
by the FOSEP approximation but not by the RPA. In the present simple model this goes even beyond second order 
as can be seen from Eq. (J4q). 



-1 



>> 
'a 

x 

CO 



triplet 






FOSEp(p) 




" TDA 


0.2 


0.4 0.6 0.8 


singlet 


FOSEp(p) ^^-"""^ 
^^^^ TDA - 






FOSEP(U) 



0.2 0.4 0.6 0.8 1 
perturbation parameterA 

FIG. 1. The eigenvalues of the FOSEP, RPA, and TDA equations for the triplet and the singlet excitation as described in the 
text. The physical FOSEP excitation energies AE FOSBP( - p) (labelled as fosep(p)) coincide with the exact excitation energies. 
The parameters of the Hubbard model are chosen as t = 1 and U = 10. All energies are given in units of £. For the triplet 
excitation, the RPA becomes instable yielding complex eigenvalues for A > 0.2. 

For the singlet excitation, results analogous to the triplet case are found. Again the FOSEP result AE S " K ' 
coincides with the exact excitation energy AEs u — Es u — E . The TDA result evaluates to 



AE s° A 



= 2t+ X -U 



(51) 



while the RPA yields 
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ae ^Mp) = ^4 t 2 + 2 \tU. (52) 

A plot of these result with the same choice of parameters as for the triplet excitation is shown at the bottom of Fig. 1. 

The Hubbard parameters used for the plots in Fig. 1 have been chosen such that the differences between the three first 
order approximations are large and become apparent. When the on-site repulsion term U is decreased in comparison to 
the hoppingparameter t, the RPA becomes stable and the differences between the discussed approximations diminish. 
From Eq. ([44|) it can be seen that U is the source of ground state correlation in the model. Thus the present analysis 
supports the findings of Sec. II D| and leads to the conclusion that among the three considered approximations for 



excitation energies only the FOSEP method includes ground state correlation in a consistent manner. 

Of course, it is a special peculiarity of the present simple model that the FOSEP approximation already yields 
the exact excitation energies. Certainly this model may seem inadequate to draw conclusions on the performance of 
FOSEP for realistic large finite quantum systems which constitute the main area of possible applications. Therefore, 
this model should not be understood as a test for large-scale numerical calculations but rather as the simplest possible 
system where the effect of ground state correlation could be investigated analytically. 



IV. CONCLUSIONS 



In this paper we have presented a novel approximation scheme for the calculation of energies and transition moments 
of many-body systems. It derives from the first order approximation to the self energy of a recently proposed extended 
particle-hole Green's function fl4|| . The resulting approximations for excitation energies and transition moments prove 
consistent in first order but also higher order terms are present. Starting out from a Hartree-Fock single-particle 
description, the FOSEP approximation yields a matrix eigenvalue problem of the same size as the well known Random 
Phase Approximation (RPA). In contrast to the RPA, however, FOSEP presents a hermitian eigenvalue problem and 
thus avoids the instabilities of the RPA. Although the FOSEP approximation has many properties in common with 
the RPA, like size consistency and the invariance with respect to unitary transformations of the Hartree-Fock virtual 
orbitals, there are also substantial differences. We have shown by a perturbation theoretical analysis up to second order 
that the FOSEP approximation for the excitation energies consistently includes part of the ground state correlation 
whereas the RPA proves inconsistent in this respect. This statement is supported by the results of a very simple 
exactly solvable model. In the considered model the ground state is correlated whereas the approximated excited 
states are not. It turned out that the FOSEP approximation yields the exact results for the excitation energies 
whereas the RPA does worse than the simpler Tamm-Dancoff approximation (TDA). We also have addressed the 
equivalence of length and velocity form of the transition moments. The exact preservation of this equivalence is a 
very peculiar property of the RPA. Within the FOSEP approximation this equivalence is only preserved in first order 
which is consistent for a first order scheme. The second order error encountered, however, may provide useful in 
estimating the applicability of the approximation. 

Future calculations on realistic systems still have to provide the ultimate test for the usefulness of the approximation 
scheme presented here. From the present point of investigation the FOSEP method seems to have excellent prospects 
to find widespread application like, e. g., in clarifying the electronic structure of larger molecules especially when 
ground state correlation is important. A possible line of extending the FOSEP method is to start out from a multi 
configurational self consistent field (MC-SCF) reference state instead of the Hartree-Fock Slater determinant \Q ). 
This would allow for an adequate treatment of open-shell or other systems where strong ground state correlation 
prohibits the zeroth order description by a Slater determinant. The first order particle-hole self energy seems ideally 
suited for such an extension since the primary matrices are defined without reference to occupied or virtual Hartree- 
Fock orbitals in contrast to the RPA or TDA. Thus the inclusion of multi-configurational reference states provides a 
natural extension of the present theory. Straightforward approximations to the static particle-hole self energy can also 
be obtained from a ground state description by density functional theory (DFT) . While DFT has been very successful 
in predicting ground state properties, the DFT description of excited states is a vivid and still open field of current 
interest p| , p2| , ^3|JT^ ] . The static particle-hole self energy seems well suited for adaption to DFT since it provides a 
simple model for excitation properties that only requires a decent approximate description of one and two particle 
densities in the ground state as input as discussed in the appendix. A direct formulation of the particle-hole self 
energy in terms of density functionals, on the other hand, may as well lead to powerful approximations. Another open 
point is the development of higher order approximations to the particle-hole self energy. This will allow to increase 
the accuracy and lift the present restriction to particle hole type excitations. A realisation of systematic higher order 
approximations could follow the concepts of correlated excited states and intermediate state representations developed 
in Refs. jl^Jn]] and is left for future work. 
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APPENDIX 



The general formalism behind the theory of extended two-particle Green's functions is described thoroughly in 
Ref. [[14| where also the definition of the states |i^- s ), the metric p,, and the generalised excitation energy operator Ti. 
can be found. Here we only want to show one result of the general theory: The general expression for the primary 
block 7£ of the excitation energy operator matrix is given by 



n 



= {Y ra \fM\Y r , 



a\a r , 



H, a J r ,a s > J |* ) 
{[#,4,] , a r } \* )(* \a*al\V ) 



+ (*„! {a r ,al,} |* ><* | 



H , a s 



at |*o) 



H , al 



a r |*o)(*o| {a s > , a\} 



+<* |a r 4,|* )<* | {[.ff,<v] , at} |* ) 



(Al) 



Here |^ ) denotes the exact ground state of the system. In general, Eq. (Al) requires the evaluation of ground state 
expectation values of one and two-particle operators. This is due to the particular combination of commutators and 
anticommutators and to the fact that the Hamiltonian H is a two-particle operator. But this means that the primary 
block iL and therefore the static particle-hole self energy can be calculated exactly if the general one and two-particle 
densities of the ground state are known. Approximating the exact ground state |\I/ ) by a Slater determinant leads 
to a factorisation of the two-particle densities and the first order expression ([l2]) is obtained. Other approximations 
for the densities than those obtained by the Hartree-Hock Slater determinant are of course also possible. Density 
functional theory, on the one hand, or a multi-configurational self consistent field (MC-SCF) approximation for the 
ground state wavefunction, on the other hand, provide interesting alternative approaches. 
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